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ABSTRACT 


We describe a new atomic and molecular database we developed for use in 
the spectral synthesis code Cloudy. The design of Stout is driven by the data 
needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic 
temperatures 2.8 K < T < 10^'^ K and densities spanning the low to high-density 
limits. The radiation held between photon energies 10“® Ry and 100 MeV is 
considered, along with all atoms and ions of the lightest 30 elements, and ~ 10^ 
molecules. For ease of maintenance, the data are stored in a format as close as 
possible to the original data sources. Few data sources include the full range of 
data we need. We describe how we hll in the gaps in the data or extrapolate 
rates beyond their tabulated range. We tabulate data sources both for the atomic 
spectroscopic parameters and for collision data for the next release of Cloudy. 
This is not intended as a review of the current status of atomic data, but rather 
a description of the features of the database which we will build upon. 

Subject headings: atomic data; atomic processes; molecular data; database 



3 


Introduction 


Cloudy is an openly available spectral simulation code based on detailed microphysics, 


most recently reviewed by Ferland et al. (2013a). It considers microphysical processes from 


hrst principles to determine the excitation, ionization, and thermal properties of a mix of 


gas and dust. Much of this physics is described in Osterbrock & Ferland (2006), hereafter 


AGN3. A very wide range of densities and temperatures can be modeled, and the full 
radiation between 10“® Ry and 100 MeV is considered. 

Massive amounts of atomic and molecular data are needed to do such simulations. 
These include energy levels; transition probabilities; collision rates with electrons, 
protons, and atoms; photoionization cross sections; collisional ionization rate coefficients; 
recombination rate coefficients, along with charge exchange ionization/recombination data. 


There are several spectral databases available, including Chianti (Dere et al. 1997 Landi 


et al. 2012) LAMDA (Schoier et al. 2005), JPL (Pickett et al. 1998), and CDMS (Muller 


et al. 2001| 2005). These provide energy levels, and transition probabilities. Chianti and 
LAMDA also include collision rates with a particular emphasis on certain applications. 
Chianti and LAMDA are included in the Cloudy distribution and are used in our simulations 


Ferland et al. 

2013a 

Lykins et al. 

2013 


At times during the development of Cloudy, we have need to create additional models 
of atoms or molecules. What format should we use? Chianti comes closest to providing 
the data we need, but its format does not allow for more than 999 levels and the collision 
rates are presented in a format that is far removed from the original published form. Only 
a few spline points are given for collision rates, and they emphasize temperatures higher 
than those found in photoionization equilibrium, so the fits are sometimes not valid at the 


low temperatures we need (Ercolano et al. 2008). Chianti’s use of spline interpolation can 


lead to unphysical negative collision strengths in Chianti version 7. Furthermore, we must 
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include collisions with atoms and molecules, which are important for photodissociation 
regions (PDR) calculations. Hence the need for our own database. 

This paper describes how we implemented our spectral line database. It is not intended 
as a dehnitive reference for the state of the art in atomic and molecular data today. 
Continuous updates to the database will occur and be described in future papers. 


2. The Stout Database 

The new database was designed to have the following properties: 

• The data format must be easy for a person to maintain since continual updating is 
necessary. 

• It must provide for different types of data. For example, radiative rates might be 
specihed as oscillator strengths, transition probabilities, or line strengths. 

• Collision data should, if possible, cover the temperature range considered by Cloudy, 
currently 2.8 K to 10^° K. This is seldom available so we need to have a strategy to 
extrapolate beyond the limits of the tabulated data. 

• We must be able to reliably interpolate upon tables of rate coefficients without 
producing unphysical negative values, which may introduce negative collision rates. 

• Both resonance and subordinate lines must be included since Cloudy is applied to 
dense environments where subordinate lines are important. 

• Both molecular and atomic data must be considered. 

• A broad range of collision partners, including electrons, H 2 , H°, He, and H’*', must be 
considered. 
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• Each file must explain its provenance by documentation at the end of the hie. 

• As far as possible, the data must be presented in their original format. We use 
the tabulated collision rates, collision strengths, energies, etc, as they appear in the 
original publication. This makes the data much easier to maintain. 

• Numbers within data hies are free format. Each number need only be surrounded by 
a space or tab character to distinguish separate entries. This makes it both easier to 
maintain and to remain close to the format of the original data source. 

• There must be no limit to the number of levels in a model. 


2.1. Spectral models in Cloudy 


We begin with a description of the atomic models in Cloudy. Cloudy has two distinctly 
diherent types of atomic models due to the diherent level structures of various isoelectronic 


sequences (Ferland et ah 2013a). The H- and He-like isoelectronic sequences have excited 
states that are closer to the continuum than the ground state. As a result, these states 
are strongly coupled to the continuum, with levels populated following recombination or 
collisions between excited states or the continuum. Excited states are relatively weakly 
coupling to the ground state. The H- and He-like isoelectronic sequences use the models 


described by the series of papers starting with Bauman et ah (2005) and Porter et ah 


(2005). Porter et al. (2009, 2012, 2013) give the most recent updates. 


For more complex ions, the lowest excited states are close to the ground state and 
are strongly coupled to it. The influence of the continuum is weak. The remaining atoms 
and ions are treated with the atomic models described in this paper or with Chianti. The 
Appendix contains Table that summarizes the data sources we use. We describe how we 
use these data in the following subsections. 
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The data for the atomic species, contained in the Stont database, are located in 
separate directories for each species with a structnre similar to that of Chianti. Each data 
set consists of three hies - energy levels (the hie with the extension “nrg”), transition 
probabilities (extension “tp”), and collision data (extension “coll”). 


2.2. Energy levels 


We use the experimental level energies from NIST (Kramida et ah 2014) if possible. 
Experimental data are available for most species. The level energies given in the NIST 
database are usually derived from measured line wavelengths. If there are no experimental 
data, we utilize theoretical data which often are of lower accuracy. 


Transition rates that come from theoretical calculations must be corrected for any 
diherences between experimental and theoretical energies. The level ordering may not 
agree so it is absolutely important to match the level assignments given in diherent data 
sources. When this issue is overcome, the integrity of the particular system is assured, 
and the calculated collisional parameters and radiative parameters are consistent with the 
experimental energy levels. 


By default we report wavelengths, in Angstrom that are derived from the stored energy 
levels (so-called Ritz wavelengths). We do allow the wavelength to be specihed to override 
this default. The convention in atomic physics is to use air wavelengths for A > 2000 A 
and vacuum wavelengths for A < 2000 A. More recently, work has started to appear which 
uses vacuum for all wavelengths. The Sloan project (see The Sloan Digital Sky Survey at 
http://www.sdss.org) is an example. By default we follow the atomic physics convention 
but provide an option to report only vacuum wavelengths. The index of refraction of air is 


taken from Peck & Reeder (1972). 








7 


Table 1: A sample of the level energies file (s_2.nrg) for S II from Stout. The first column 
represents the level index that is used in the transition probability and collision data files, the 
second column gives level energies (in cm“^), the third column gives level statistical weights 
g, and the last two columns give a level designation with standing for odd-parity levels. 


N 

energy 

9 

conhguration 

LSti 

1 

0.00 

4 

3s2.3p3 

4S* 

2 

14852.94 

4 

3s2.3p3 

2D* 

3 

14884.73 

6 

3s2.3p3 

2D* 

4 

24524.83 

2 

3s2.3p3 

2P* 

5 

24571.54 

4 

3s2.3p3 

2P* 

6 

79395.39 

6 

3s.3p4 

4P 

7 

79756.83 

4 

3s.3p4 

4P 


We present a sample of a hie with columns description for the energy level data in 
Stout in Table Just a small part of the hie s_2.nrg with the S II ion level energies are 
given here. The complete data table is given by Kisielius et ah (2014a|, whereas level data 


are taken from the NIST database (Kramida et al. 2014). Only the level numbers, energies, 
and the statistical weights are utilized for deriving of the transition probability or collision 
strengths, whereas the conhguration and LSti are given for information purposes. 


2.3. Radiative transitions 

The radiative transition between the upper level u and the lower level I can be 
parameterized as a line strength S, oscillator strength / (or gf), or transition probability 
A, although only the latter enters in a calculation of level populations and emission spectra. 
Diherent data sources will provide diherent parameters, and we accept all three. 

We prefer to utilize the transition line strength S over the weighted oscillator strengths 
gf or transition probabilities A. The advantage to S is that it does not depend explicitly 
on the transition energy AE (or the transition wavelength A), whereas gf and A do. Many 
published transition data are the result of theoretical calculations and use theoretical 















Table 2: Conversion factors and coefficients from the transition line strengths S (in a.u.) to 
the radiative transition probabilities A (in s“^). 


type 

factor Cx 

C-A 

Cae 

El 

GAiT^e'^a'^ 

2.02613 X 10^® 

2.14200 X 10^° 

3h 

Ml 

647r‘^/A3 

3h 

2.69735 X 10^3 

2.85161 X 10^ 

E2 


1.11995 X 10^® 

5.70322 X 10^ 

Ibh 

M2 

647r®/AgaQ 

15h 

1.49097 X 10^3 

7.59260 X 10-1 

E3 

20487r®e2a^ 

3.14441 X 10^^ 

7.71311 X 10-2 

4725/1 

M3 

20487r*pgag 

4.18610 X 10^2 

1.02683 X 10-® 

4725/1 


energies while we use experimental energies where possible. Therefore, a correction due to 
the uncertainty in the calculated transition energy AE, or wavelength A, values must be 
done. The conversion to the experimental transition energies AT^exp or observed transition 
wavelengths Aexp is: 

= /l5.(A.k/A„p)“+‘ (1) 

where k is the transition multipole order {k = 1, 2, 3,...), A is the transition probability, Ath 
is the theoretical transition wavelength. 

The transition line strength S is expressed in atomic units (a.u.). It is symmetric in 
relation to the initial and hnal states, and is obtained as a square of the corresponding 
El, Ml, E2, M2, E?>, M3 transition matrix elements. In this case, the electric multipole 
emission transition probability (Einstein A-coefficient) A^^ (in s“^) can be determined as 

At, = C'aS/(9„A“+') (2) 

where Qu is the statistical weight of the upper level, and C\ is the conversion factor. The 
expressions of the factor C\ for various multipole transitions are presented in Table We 
provide numerical values of the conversion coefficients when A is expressed in A and the line 
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Table 3: A sample of the transition probabilities file (s_2.tp) for S II from Stont. The hrst 
colnmn represents radiative transition data type (“A” for a transition probability, “f” for 
weighted oscillator strength, “S” for a line strength), the second colnmn gives the lower level 
index, the third colnmn gives the upper level index, the fourth column gives a transition 
parameter value, and the hnal column points to the radiative transition type {El, E2, E3, 
..., Ml, M2, ...). 


data 

type 

Ni 

Nu 

TP 

transition 

type 

S 

1 

2 

5.54E-03 

E2 

S 

1 

2 

1.77E-05 

Ml 

s 

1 

3 

1.29E-02 

E2 

s 

1 

3 

6.30E-07 

Ml 

s 

1 

4 

2.25E-06 

E2 

s 

1 

4 

3.37E-04 

Ml 

s 

1 

5 

3.00E-10 

E2 

s 

1 

5 

1.67E-03 

Ml 

s 

1 

6 

2.65E-01 

El 


strength S is calculated in a.u. In the case when we have the transition energy AE instead 
of wavelength, we can use similar expression for 

At, = CiE(Ai5)“+'S/9„. (3) 

The conversion coefficients Cae are given in Table for AE determined in a.u. For 
different energy units, one must rely on these standard relations: 1 a.u.= 2 Ry= 27.211385 
eV= 219474.63 cm“\ and the inverse relations: 1 Ry= 0.5 a.u.; 1 eV= 0.036749324 a.u.; 

1 cm“^ = 4.5563353 ■ 10“® a.u. Having the radiative transition probabilities of various 
multipole orders, one can simply derive the absorption oscillator strengths fiu using a 
simple expression: 

flu = 1.4992 X 10-^^\‘^{gu/gi)Aui (4) 

where gi is the statistical weight of the lower level. The relation between the oscillator 
strength / and transition probability A does not depend on the radiative transition type. 

Table gives a sample of the transition data file for S II, with the data coming from 
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Kisielius et al. (2014a). The transition line strengths S are given as the basic radiative 


transition data as they do not depend explicitly on the transition energy. The NIST 
database traditionally provides the radiative transition probabilities (rates) A. Conversion 
from the transition line strengths S to the transition probabilities A depends on the 
transition type. It can be performed with the help of Table 


2.4. Collisions 

Collisional data can be given as collision strengths, effective collision strengths, collision 
cross sections, and rate coefficients. In the electron-impact excitation: 

A^^{Ei) + e{£i) —t A^~^[Eu) + e(£u), (5) 

the energy conservation law leads to 

El + £i = Eu + £u- ( 6 ) 

Here Ei and Eu are the energies of the lower and npper levels, Si and Su are the kinetic 
energies of the incident and the scattered electron. 

Maxwellian-averaged effective collision strengths are ntilized in Stont. Here we provide 
the basic relations for these parameters, whereas data sonrces are given in Table Onr 
preferred method is to use collision data calculated in the close-coupling (CC) approach, 
e.g., the R-matrix method. Unfortunately, such data are not available for many ions. Even 
in the case when some data exist, these usually deal with only a few of the lowest levels 
or even the LS terms with unresolved fine-structure levels. So one must resort to less 
elaborate approaches, such as the distorted-wave method, the plane-wave approximation, 
or a ^ (g-bar) formula, for the remaining data. 

The dimensionless collision strength hi is the best to describe the electron-impact 
excitation process from the lower level Ei to the excited level Eu- It is symmetrical in 
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regard to the initial and final states parameter, i.e, Qiu = ^ui- For ions, it has a finite value 
at the excitation threshold and varies only slightly with the incident electron energy if 
autoionization resonances are not considered. For neutral atoms, the collision strength Qiu 
goes to zero at the excitation threshold. 


At high incident electron energies, the behavior of the collision strength depends on 
the transition (line) type and has a different form for allowed or forbidden transitions. The 
collision strength fliu is determined as a square of the excitation operator’s matrix element. 
It is connected to the excitation cross-section aiu and the de-excitation cross-section aui by 


simple relations: 

aiu{ei) = (7) 

giEi 

and 

aui{e^) = ( 8 ) 

Qu^u 

where ttuq = 8.7972 x 10“^^ cm^ is the atomic cross-section unit. For electric dipole allowed 
transitions, one can express the excitation cross section aiu through an effective Gaunt 


factor giui^i) (as in Mewe 1972), often called the “g-bar approximation”: 

/ \ flugi') 

- V3e,(E. - E.) 

with flu being the absorption oscillator strength. 


( 9 ) 


The collision strengths are integrated over a Maxwellian distribution of free-electron 
energies in order to determine the effective collision strengths Tiu = Tui (or rate parameters) 
at some electron temperature Tg: 

poo 

'^iu= D.iui£u)exp{-eu/kBTe)d{-eu/kBTe). ( 10 ) 

Jo 

Here /cb refers to Boltzmann’s constant. In this case, the excitation rate coefficient qiu (in 
cm3 jg expressed as 

gexp((F;; - Eu)/kBTe)Tiu 


qiu = 8.629 X 10 


( 11 ) 
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Table 4: A fragment of the effective collision strength hie s_2. coll for the ion S II from 
Stout. The hrst column represents data type, e.g., “TEMP” stands for a temperature grid 
(in K), “CSELECTRON” for the effective electron-impact excitation strength T, “RATE 
PROTON” for the proton excitation rate, etc. The data types and their sources are provided 
in the same hie in the comments lines. The second column gives the lower level index, 
the third column gives the upper level index, and next columns give a particular collision 
parameter for the corresponding temperature. 


data 

type 

Ni 

Nu 






TEMP 



5.00E+03 

7.00E+03 

I.OOE+04 

I.50E+04 

2.00E+04 

CSELECTRON 

1 

2 

2.66E+00 

2.62E+00 

2.56E+00 

2.48E+00 

2.4IE+00 

CSELECTRON 

I 

3 

3.98E+00 

3.9IE+00 

3.83E+00 

3.7IE+00 

3.6IE+00 

CSELECTRON 

I 

4 

6.86E-0I 

6.94E-0I 

7.04E-0I 

7.I7E-0I 

7.27E-0I 

CSELECTRON 

I 

5 

I.38E+00 

I.39E+00 

I.42E+00 

I.44E+00 

I.46E+00 

CSELECTRON 

I 

6 

2.25E+00 

2.36E+00 

2.54E+00 

2.75E+00 

2.84E+00 

CSELECTRON 

I 

7 

2.0IE+00 

2.09E+00 

2.I9E+00 

2.30E+00 

2.3IE+00 


Note. — Either deexcitation rate coefficients (cm^ s or effective collision strengths can be specified. 
The colliders include electrons, protons, alpha particles, He’*', He°, H2 (ortho and para), and H°. 


whereas the de-excitation rate coefficient qiu is determined by formula: 

Tiu 


qui = 8.629 X 10’ 


SuTe 


1/2- 


( 12 ) 


Table 1^ gives a fragment of the collision data hie s_2.coll for S II. For this ion we employ 


the electron-impact excitation data from Tayal & Zatsarinny (2010). There can be several 


temperature grids in one data hie, especially when diherent projectiles, such as electrons, 
protons, hydrogen atoms, or hydrogen molecules, are described. In a similar way, there can 
be diherent data sources for diherent level combinations even for the same collider. The 
collision data for a given transition and collider will be overwritten if new data for that 
transition and collider appear later in the hie when read in by Cloudy. 


Most of the data sources in Table are given for the electron collisions. Nevertheless, 


some colliders other than electrons are included. For an atomic hydrogen collider, very 
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important in PDRs, we use data from Launay & Roueff (1977) for C I, from Abrahamsson 


et al.l (2007); Krems et al. (2006) for O I, from Hollenbach & McKee (1989) for Ne II, 


from Barklem et al. (2012) for Mg I, and from Barinovs et al. (2005) for Si II. For proton 


colliders, very important in collisions between levels with similar energies, we utilize data 


from Roueff & Le Bourlot (1990) for C I and data from Pequignot (1990) for O I. For H 2 , 


data from Schroder et al. (1991) are employed for C I, data from Wiesenfeld & Goldsmith 


(2013) for the C II ion, and data from Jaquet et al. (1992) for O I. For helium impact on 


neutrals, we use data from Staemmler fc Flower] (1991) for C I, and from Monteiro & Flower 


(1987) for O I. 


For collisions involving molecules, the literature often gives deexcitation rates rather 
than collision strengths. We accept deexcitation rates for any transition and species. 

It is necessary to interpolate within tables of collision rates versus temperature, and 
in many calculations, extrapolate beyond the tabulated range. Within the table, we 


interpolate using the method of Fritsch & Rutland (1984), which is local and piecewise 


cubic, and maintains the monotonicity properties of the underlying data. This ensures that 
the interpolation does not introduce any “overshoots”, where the interpolated value does 
not he within the range of the tabulated data. Such overshoots appear to be the source of 
the negative collision strengths that are present in version 7 of the Chianti database. 


2.4-1- Temperature extrapolation for atoms and ions 

Cloudy considers the temperature range extending from 2.8 K to 10^° K and considers 
all ions of the first thirty elements along with several dozen molecules. Gaps in the collision 
data are common. Often we must extrapolate beyond the range of the tabulated data, or 
improvise entire collections of data. 

For temperatures below the range of the tabulated data, and for ions with positive 
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charge, we assume that the collision strength is constant to extrapolate below the lowest 
tabulated temperature. Physically, an effective collision strength is a Boltzmann average 
over the excitation cross section. As the temperature goes to zero, this average is over a 
narrow range near threshold, and will tend to be constant. This is not true if there are 
strong resonances very near threshold but it is a reasonable first approximation. For neutral 
species the collision strength goes to zero at energies near threshold so the effective collision 
strength also goes to zero as T —)■ 0. We do a linear interpolation between the lowest 
temperature value and 0.0. We use these collision strength laws to form the appropriate 
temperature scaling when working with data giving collisional deexcitation rates. 


For high temperatures we use Burgess & Tully (1992) to guide the extrapolation. 


Burgess & Tully (1992) consider three possible types of transitions with different behavior 


at high energies (temperatures). Type 1 for the electric dipole transitions. Type 2 for the 
non-electric dipole, non-exchange transitions. Type 3 for the exchange transitions, with 
Types 4 and 5 being special cases (for more details on the transition classification see 


Burgess fc Tully (1992)). We use the first two types. Type 1 and Type 2. The Type 3 


classification can be avoided when levels and transitions are expressed in the intermediate 
coupling rather than pure LS coupling. In this case one can not separate the spin-changing 
transitions since the selection rules are applied for the total angular momentum J. 


Transition types can be deduced from the energy levels files (*.nrg) where J, 
configurations and their parities are given or from the radiative transition files (*.tp) where 
transition types are given (but these are not present in all transition data files). In general, 
our data files contain information necessary to make a separation between Type 1 and Type 
2 transitions. For the Type 1 transitions, the high-temperature behavior of the effective 
collision strength T is described by a simple relation T = Ci In(Te). The value of Ci can 
be derived from the last tabulated temperature point in the data file. For the Type 2 
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transitions, the effective collision strength does not depend on the electron temperature Tg, 
i.e., T = C 2 . The value of C 2 is the value of T at the last tabulated temperature. 


Tests show that the low-temperature extrapolation does affect calculations. In 


photoionization equilibrium very low kinetic temperatures are possible (Osterbrock & 


Ferland 2006). The constant temperature test cases in the Leiden PDR comparison (Rollig 


et al. 2007) have Tkm = 50 K, lower than many tabulated rates. Predictions of some Leiden 


test cases were affected by the form of the low-T extrapolation. 


2.4-2. Gaps in the collision data 


We use the g approximation (Seaton 1962 van Regemorter 1962) to fill in missing 


electron collision data. This is a highly approximate relationship between the transition 


probability or oscillator strength, and the collision strength. We use Mewe (1972) for those 


isoelectronic sequences he considered, and van Regemorter (1962) for others. 


We provide a way to test the effects of such uncertain data. Cloudy includes a built-in 
Gaussian random noise generator. This was used, for instance, to assay the effects of 
missing H 2 collision rates upon the final spectrum ( Shaw et al.||2005 ). Repeated calculations 
will reveal the uncertainties introduced by the approximations, if the uncertainties can be 
quantified. 


Some databases have no radiative transition between large blocks of levels. For 
instance, a species may have no El transitions between the ground and first excited 
configurations. Higher order transitions are possible but many databases present only El 
transitions. If theoretical collisional rates have not been computed, then there would not 
be any coupling between the configurations. It is not possible to simultaneously solve for 
the populations; the matrix becomes ill conditioned. In cases where we have no radiative 
or collision data, we leave the radiative transition rate as zero and use an electron effective 
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collision strength of 10 . This was chosen to be as small as possible while allowing the 

linear algebra to function properly. 

2.4.3. Temperature extrapolation for molecular excitation 

When molecular collisional deexcitation rate coefficients q{T) are provided only over a 
limited temperature range, the following two simple extrapolation approaches are applied: 

g(T) = g(Tiow), T < Tiow, (13) 

and 

/ T \ 

g(T) = g(Thigh) exp(-T/10' K), T > Thigh- (14) 

V^high/ 

Tiow and Thigh correspond to the low and high temperature limit, respectively, of the data. 
The extrapolation formulae are valid for inelastic collisions of neutral molecules (e.g., OH) 
or molecular ions (e.g., HCO’*') with neutral colliders (e.g., H, He, H 2 ) for deexcitation 
(downward) transitions resulting in changes in hue-structure, rotational, and/or vibrational 
levels of the target molecules. Physical justihcations for the extrapolations as well as 
caveats for their use are described below. 

Extrapolation to low temperature 

The deexcitation rate coefficients as a function of temperature T are obtained by thermally 
averaging the inelastic integral cross sections over a Maxwellian kinetic energy distribution 
given by 

( 8 k T\ 1 r'^ 

- ~) tl / 0-n^z(^kin)eXp(-Tkin//CB7')TkindTkin, (15) 

Trp J [k^Ty Jo 

where (j„_s.«(Tkin) is the state-to-state inelastic cross section, Tkin the center of mass kinetic 
energy, p the reduced mass of the collision complex, and I {u) the lower (upper) levels in 
the molecule. 
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Rewriting Eq. (15) with the cross section in terms of the relative velocity v of the 


collision system gives (Flower 1990) 


1/2 ^ ^ 3/2 ^oo 

TT 


/ JO 

If the cross section is assumed to have the analytical form 

au^i{v) = 


au^i{v) ex.p{—iJv‘^/2kBT)v^dv. 


(16) 


(17) 


for all V (or E kin ) where B is an (undetermined) constant and a is some power, then the 
rate coefficient takes the form 

g(T) = A{a)B{T/fif (18) 


(Walker et ah 2014). Here b = {a + l)/2 and H is a function of a, both deduced from the 
Gaussian integral in Eq. (16). This result is exact, given the assumption of Eq. 0 , and 
applicable to all collision systems. It is approximate if the cross section dependence varies 
with V as in real systems. 


Now at sufficiently low kinetic energy, Wigner (1948) showed that the inelastic cross 
section takes the form 

~ dr/''/ (19) 

where £ is the total orbital angular momentum of the collision complex. In most systems 
of interest, s-wave scattering (i.e., £ = 0) is allowed and dominates at low kinetic energy. 
Therefore, 

( 20 ) 


O'u^l ~ EyA ~ V 


-V2 _ , ^,-1 

kin 


or a = —1, & = 0 and the rate coefficient becomes a constant, independent of temperature 


as given by Eq. (13). 


Equation (13) is absolutely valid under two conditions: i) when T < Twigner, the 


so-called Wigner regime, where higher partial waves {£ > 0) do not contribute to the cross 
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section and ii) when relevant selection rules do not forbid s-wave scattering. Typically, 


Tw\ 


Wigner 


niK and much less than Tiow However, for T ~ 1 — 100 K the rate coefficient 


is usually oscillatory due to the presence of orbiting and Feshbach resonances in the cross 
section. This behavior cannot be easily analytically reproduced so that extrapolating the 


Wigner threshold behavior of Eq. (20) to Tiow is a reasonable pragmatic approach. 


In the event that s-wave scattering is forbidden for the particular transition (which is 
rare), the cross section would drop rapidly to zero as v goes to zero 


~ ~ V. 


( 21 ) 


The rate coefficient would be overestimated by Eq. (13), but this error would be limited to 
low astrophysical temperatures, T ^ 10 K. 


Extrapolation to High Temperature 


A number of approaches have been proposed for extrapolating the deexcitation rate 


coefficient to higher temperatures beyond Thigh- For linear molecules, Schoier et ah (2005) 
htted the available data in LAMDA to the form 


9(T) = ^exp(-,3/(«:BT)‘''*)exp(-7/(fer)‘/=) 


( 22 ) 


where a, /I, and 7 are £t parameters and then used the £t for T > Thigh- A more pragmatic 
approach, which avoids htting, is to apply a hard-sphere model- This assumes that the 
cross section is independent of kinetic energy giving a = 0 and & = 1 in Eqs- ( p!^ and (18), 
respectively, so that 

q{T) ~ T^A. (23) 


However, the inelastic cross section typically turns up to a maximum near a few eV before 
decaying at higher energies due to the increasing importance of collisional dissociation, 
electronic excitation, and collisional ionization. Therefore, to prevent the rate coefficient 
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from growing too large at high T, the relation (23) is multiplied by an exponential damping 


factor to give Eq. (14). The exact form is not important as the molecular abundances 
decrease rapidly for T > 5000 K. 


3. Other details 

3.1. Baseline models, unmodelled species 

Many species have level energies and transition probabilities tabulated in NIST, but 
have no electron collisional rates at all. For these species we created “baseline” Stout data 
files. These contain the NIST level energies and transition probabilities but use the g 
approximation for all collision data. These are marked as “baseline” in Table 

It was not possible to create models for all ions. NIST did not have sufficient data to 
compute models for the following species: F I, Cl XIII, Cl XV, Sc I, Ti I, Ti II, V I, V II, 
V III, V V, V XIV, Cr I, Cr III, Cr V, Mn II, Mn III, Mn IV, Mn VII, Co I, Co IV, Co V, 
Co VI, Co VII, Co IX, Ni VI, Ni VIII, Ni X, Cu II, Cu III, Cu IV, Cu V, Cu VI, Cu VII, 
Cu VIII, Cu IX, Cu X, Cu XII, Cu XII, Cu XIX, Cu XX, Zn I, Zn III, from Zn V to Zn XIX, 
and Zn XXI, Zn XXII, Zn XXVI. As a result Cloudy calculations do not predict lines of 
these ions. Calculating sufficient data for these species should be a high priority. 

3.2. Masterlists - specifying which database 

Cloudy uses a total of three databases. Stout, described here, along with Chianti and 
LAMDA. Each database has its own “masterlist”, a file specifying the species present in 
that version of the database. The Stout masterlist file was used to derive Table HI 

In a particular Cloudy calculation, each of these masterlist files will be read. It is 
likely that a particular species is present in more than one database and its masterlist file. 
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The priorities for deciding which database to use are: 1) the H- and He-like isoelectronic 
sequences are always treated with our unihed model, 2) Stout, 3) Chianti, and 4) LAMDA. 


3.3. Suprathermal electrons 

When cosmic rays or Auger electrons enter neutral gas they create a population of 
suprathermal electrons which ionize and excite the gas ( Spitzer fc Tomasko||1968 ). We solve 
for the population of these suprathermal electrons explicitly. 


These electrons have an energy of typically 20 - 40 eV and can cause internal excitations 
of all atoms and molecules. We include this as a general excitation process using the Born 


approximation outlined by Shemansky et ah (1985) 


4. Discussion and Summary 


As described in the mandate for the development of Stout, it is now far easier to 
maintain and update the line database in Cloudy, and to add entirely new species. With 
the addition of the species given in Table there are now far more lines predicted than in 
previous versions of the code, producing far richer spectra. 

Figure shows an example. This is a coronal equilibrium metal-rich (Z = 35Zq) gas 
with a density of 1 cm“^ and a temperature of 5 x 10^ K. The panels compare the current 
version, soon to be released as C15, with CIO, the last Cloudy release before beginning the 
move of ionic models to external databases. There are now a far greater number of faint 
lines. 


Despite the large increase in the number of lines the total cooling of the gas is relatively 


unaffected. Lykins et ah (2013) describe our calculation of the total gas cooling along with 


our strategy for determining how many levels to model. The cooling is dominated by a 
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few very strong lines so the large number of faint lines do not increase it significantly. The 
faint lines can be important when abundances are non-solar (as in Figure [^, low-abundant 
species are of interest, or if a number of faint lines blend to produce a stronger feature. 

This paper is a dehnition of our database and explains how Cloudy uses it. It is not 
intended as a review of the state of the art of atomic data in 2015. Future papers will 
expand the atomic/molecular data using the framework outlined here. 

The database is designed to be easy to maintain and modify due to the need to 
constantly modify it as new data appear. The format follows the original data papers as 
closely as possible. 

The methods we developed to £11 in missing data are described. The data needs of 
Cloudy are vast due to its very wide range of applicability. We frequently encounter cases 
where collisional rates are not available at all, or we need to extrapolate beyond the range 
of computed data. The g approximation is used to provide missing electron collision data. 
This approximation has a very broad dispersion and we provide a method of checking on 
its impact on predictions. When rates or collision strengths are available, but we need 
to extrapolate beyond the range of tabulated temperatures, we use physically motivated 
asymptotic limits. Tests show that predictions are mainly affected by the form of the low-T 
extrapolation. 

The Stout database is part of the Cloudy distribution, available on www.nublado.org. 
Its version number is the same as the Cloudy version number. This paper is the defining 
documentation of Stout and should be cited if the database is used outside of Cloudy. 

We thank the referee for a very helpful review of our manuscript. GJF acknowledges 
support by NSF (1108928, 1109061, and 1412155), NASA (lO-ATPlO-0053, lO-ADAPlO- 
0073, NNX12AH73G, and ATP13-0153), and STScI (HST-AR-13245, GO-12560, 


- 22 - 



Wavelength [microns] 


Fig. 1.— The predicted spectrum of a Z = 35 Zq collisional gas with a temperature of 
5 X 10"^ K. The upper panel shows the current results, for CIS, while the lower panel is for 
CIO, the last version before the move to external databases for ions. The density of lines is 
now far greater. 

HST-GO-12309, GO-13310.002-A, and HST-AR-13914) and thanks to the Leverhulme 
Trust for support via the award of a Visiting Professorship at Queens University Belfast 
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(VPl-2012-025). RK research is funded by the European Social Fund under the Global 
Grant measure, project VP1-3.1-SMM-07-K-02-013. 


A. Data sources 

This Appendix describes the data sources currently used by the development 
version of Cloudy. Species that are not explicitly listed in this Appendix use the 
Chianti database. With the combination of these data, Chianti, and our special 
treatments of the H and He-like iso-electronic sequences. Cloudy includes 
spectral models of all ions of the lightest thirty elements. 
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Table 5. Atomic data sources in Stout 


Z Species 


Data Source 


3 Li I 

4 Bel 
Be II 

5 Bl 
Bii 
Bill 

6 Ci 


Cii 

Cm 

7 Ni 
Nv 

8 Oi 


On 

9 Fi 
Fii 
Fill 
Fiv 
Fv 
Fvi 
Fvii 

10 Nei 
Nell 

11 Nal 
Nall 

12 Mg I 

Mg III 

13 All 
Aim 

Aliv 

Alvi 

14 Sii 
Sill 

Sim 


Si IV 
Si VII 
Si IX 
15 Pi 
Pii 
Pm 
Piv 


baseline ■ 

baseline 

baseline 

baseline 

baseline 

baseline 


see text 


Launay fe KouCTf 


Berrmgton et ah 


Johnson et al. 


(|1987K 


Mendoza 


Kouelt ^ Le Bourlof 


'77^ 


119831 


ychrbder et al. 




Goldsmith et 


iTWT' 


(199011 


ytaemmler Flower 


vViesenfeld Goldsmith 






19911 


120T3 7 


T9H5 


h'lscher i'achiev r 


Jaquet et al. 


120001 


Liang Badnell 

2( 

jii 



Bellet al.l 11998 

Wang & McConkev (11992 


BarklemI ||2UU7 i, 

Abrahamsson et al. i 2UU7||, 


Krems et al. 

200 


Monteiro h lower 

00 

|l 


TW 2 i, 


Kisielius et al. 

baseline 


i pUDQ ' 


' h'ischer Sz i'achiev 


{ 20041 


Butler Sz Zeippen 

(|1994 

baseline 

Lennon & Burke 

19941 


baseline 

baseline 

baseline 


Griffin et al. 


< 20011 , 

Verner private comniumcation 
baseline 


Hollenbach & McKee 


(|l989i 


Barklem et al. 


Mendoza 

(11988 

' 

Liang & 

Badne 

1 

2010 


j2012 I, Leep Sz Gallagher (1976 I, 


Dufton Sz King 

ston 

19871 

bampson et al. 

1199 

tl 


Hollenbach &: Me 


Butler Sz Zeippen 


||l994t 

vee iffi 


Barinovs et a 


Uutton 


i 


7989' 


((M)^ 


Kingston (1994 


Dutton 6z Kingston 1 1989 ^ 


Gallegari 6z ITigueiros 1 T998 \, 


Dutton et al. 


Liang et al. 




2UU9i 


Butler Zeippen 


Lennon 6z Burke 
baseline 


19941 




Krueger Sz Czvzak 

(1970 

Krueger 6z Gzyzak 

(1970 


baseline 
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Table 5—Continued 


Z Species Data Source 


16 


17 


18 


19 


20 


21 


22 


P VI 

Si 

Sii 

Sill 
Cl I 
Civ 
Cl VI 
Civil 
Cl VIII 
Cl IX 

An 
Am 
Arm 
Ar IV 
Arv 
Ar VI 
Ki 
Kii 
Km 
Kiv 
Kvii 
K vm 
Kx 
Cai 
Cam 
Caiv 
Ca VI 
Sc I 
Sc II 
Sc III 
Sc IV 
Sc V 
Sc VI 
Sc VII 
Sc VIII 
Sc IX 
Sc IX 
Sex 
Sc XI 
Sc XII 
Sc XIII 
Sc XIV 
Sc XV 
Sc XVI 
Sc XVII 
Sc xvm 
Tim 
Tiiv 
Ti V 
Tivi 
Tivii 


baseline 


Hollenbach & McKee 

(1989 

Kisielius et al. 

j20i4a 


laval Zatsarinny 

010 1 

Hudson et ai. 

|2012bl 


Hollenbach 

McKee 

(1989 


baseline 


Liang et al. 

(^2009 


Liang Ha 

cfnell 

2010 

Herrington et al. 

1998 


Pelan Si Berrington 

11995 

Galavis et al. 

(1935 


Kamsbottom et al. 

19971 

Galavis et al. 

11995 


yaraph Si yto 

rey p996l 


baseline 


Pelan Si Berrington 

11995 

Galavis et al. 

(1995 



yaraph Si ytorey 

1|T996 

1 


baseline 

baseline 

baseline 


Pelan & Berrington 1 1995 I 

baseline 

baseline 


Wasson et al. 

baseline 

baseline 


(20111 


Pelan & Berrington 1 1995 I 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 


Saraph & Tully 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 


(1994 \ 


Pelan & Berrington 1 1995 I 
baseline 
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Table 5—Continued 


Z Species Data Source 


23 


24 


25 


26 


27 


28 


Ti VIII 

Ti IX 

Tix 

Tixiii 

Viv 

Vvi 

Vvii 

V VIII 
Vix 
Vx 
Vxi 
Vxii 

V XIII 
Vxv 

V XVI 

V XVII 

V XVIII 

V XIX 
Vxx 

V XXI 
Cm 
Criii 
Criv 
Crv 
Crx 
Crxi 
Crxii 
Crxv 
Mni 
Mnv 
Mn vi 
Mn xi 
Mn xii 
Mn xiii 
Mn XIV 
Mn XVI 
Fel 
Fell 
Felll 
Fe VII 
Coo 
Com 
Co VIII 
Cox 
Co XI 
Co XII 
Co XIII 
Co XIV 
Co XV 
Co XVI 
Co XVII 
Nil 


baseline 

baseline 

baseline 

baseline 

baseline 

baseline _ 

Pelan & Berrington |jl995 I 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 


Berrington et al. (1998 I 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 


Grieve &: Ramsbottom ||2012 I 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 
baseline 


Hollenbach &; McKee 

19891 

Verner et al. 

119991 

l^hang (11996 

1 , Kurucz 

(^2009 

Witthoett & 

iriadnell I 

ho8i 


baseline 

baseline 

baseline _ 

Pelan & Berrington ljl995 I 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 


Hollenbach &: McKee (1989 I 
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Table 5—Continued 


Z Species 

Data Source 

Nill 

Nilll 

Niiv 

Niv 

Nivii 

Niix 

Nixvii 

Cassidy et al. (2011 \ 
baseline 

baseline 

baseline 

baseline 

baseline 

Assarwal et al. (2007 1 , 

Hudson et al. (2U12ai 

29 Cui 
Cuxiii 
Cuxiv 
Cuxv 
Cuxvi 
Cuxvii 
Cu XVIII 
Cuxxi 
Cuxxii 
Cuxxiii 
Cu XXIV 
Cuxxv 

30 Znii 

Zniv 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

baseline 

Kisielius et al. l 2014b \ 

baseline 
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